Distance dependence of entanglement generation via a bosonic heat bath 
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Within a generaUzed Caldeira-Leggett model we analyze the conditions under which a bosonic 
heat bath can entangle two microscopic quantum systems at a distance r. We find that the attainable 
entanglement is extremely distance-sensitive. Significant entanglement can only be achieved if the 
systems are within a microscopic distance that is of order of the cut-off wavelength A of the system- 
bath interaction. At larger distances the maximal entanglement is exponentially suppressed with a 
decay length of order A. We conclude that entanglement generation via a heat bath is not suitable 
for entangling remote objects. 
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Establishing and preserving quantum-mechanical en- 
tanglement [ll between two remote microscopic physical 
systems is an experimentally challenging undertaking. 
Generally speaking, the difficulties arise from the fact 
that the systems need to significantly interact with each 
other (in order to build up an entangled common state), 
and at the same time the two systems must be thor- 
oughly shielded from interaction with external degrees of 
freedom (in order to preserve the entangled state against 
decoherenceQ). In most physical situations these two re- 
quirements appear to be contradicting to each other and 
therefore can be only partially fulfilled. 

During the last years a fresh look at this entangling 
dilemma has emerged from theoretical work on the dy- 
namics of entanglement in open systems, notably from 
the work of Braun and Benatti et al. 0,01 • It has 

been shown that under suitable conditions two two-level 
systems 0, 0, S S Q or two harmonic oscillators [1, Q 
can become entangled by mere interaction with a com- 
mon bosonic heat bath, without any direct interaction 
between the microscopic systems. In such a situation 
the coupling to the heat bath has two relevant effects: 
it leads to decoherence, as it usually does, but it also 
mediates an effective interaction between the systems. 
When the latter one is strong enough to ovcrcompcnsatc 
the decohering effect, the coupling to the heat bath may 
eventually lead to entangled microscopic systems. Entan- 
glement generation via environmental modes is, needless 
to say, a sophisticated mechanism. Its theoretical analy- 
sis therefore necessarily has to rely on idealizing assump- 
tions, and it is still not clear to which extent these as- 
sumptions can be met in real systems. 

We pursued research that especially addresses the role 
of the spatial distance r between the microscopic sys- 
tems on the entangling mechanism. In doing so, we fully 
account for dissipative system-bath interactions by in- 
vestigatingan exactly solvable model along the lines of 
Ullersma [IC 
studies 0,13 



cornparatively moderate power-law dependence observed 
in [j, 0] , here we find significant entanglement between 
the systems only if the distance r does not much exceed 
the cut-off wavelength A of the system-bath interaction. 
At larger distances, E,nax, the maximum attainable log- 
arithmic negativity (as a measure of entanglement [l^). 
decreases exponentially with a decay length of order A. 
We argue that A will be typically of order of the spatial 
extension of the microscopic systems and thus conclude 
that entanglement generation via a heat bath is limited 
to truly microscopic distances only. 

In our model, the two remote microscopic quantum 
systems are represented by two identical harmonic oscil- 
lators located on a line at positions xi/2 = ±r/2. The 
choice of harmonic oscillators makes the model exactly 
solvable. Nevertheless, we still expect it to capture the 
basic physical behaviour of any system with a discrete 
spetrum (cf. The oscillators, henceforth called the 

system oscillators, have mass m and frequency a-nd 
Pi, Qi and P2, Q2 denote their canonical variables. They 
are coupled to an extended, one-dimensional heat bath 
consisting of symmetric (oc cos kx) and antisymmetric 
(oc sin kx) harmonic modes of wavenumbers fc > and 
frequencies ujk = ck, c being the velocity of sound/light. 
Let p^^", qi/"' denote their respective canonical variables. 
Wc consider a bilinear system-bath interaction Hi, where 
the two oscillators locally couple to symmetric and anti- 
symmetric bath-modes in the same manner. 



Hi = ^gkiQi + Q2)qkC0SY + 9kiQi 



Q2)ql sin y 



I, and Caldeira and Leggett Existing 
of the distance dependence are either con- 
fined to a dissipation- free spin-boson model or treat dis- 
sipation on a perturbative level only. In contrast to the 



The coupling strengths gk may be characterized as usual 

2 

by a spectral function J{lo) := '^^k 2mltjr'^(^"^'^)- Here, 
we assume J{ijj) to be linear for small oj with a Drude cut- 
off, J{u) ~ uj ^^'^^2 , leading to ohmic damping with 
a damping constant 7. Typically, the cut-off frequency 
r2 = 27rc/A is not some intrinsic frequency of the bath, 
rather, it will be determined by the physics of the system- 
bath coupling. Then, since generally \gk\ markedly de- 
clines when falls below the spatial extension I of the 
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microscopic systems, a good order of magnitude estimate 
is A ~ /, meaning that fl ~ 2itc/1. 



We also include a counter- term Vc — J2k 



9k 



2Q1Q2 cos(fcr) + Q2) in the total Hamiltonian. Its pur- 
pose is twofold: firstly, it removes the frequency renor- 
malization caused by the coupling to the bath sec- 
ondly, it ensures that the Quantum Langevin Equations 
(QLE) which we are going to derive below will only con- 
tain retarded couplings between the oscillators. 

The dynamics of the system oscillators can be ap- 
proached by means of the Heisenberg equations of mo- 
tions for their coordinates Qi, Q2- Following the analysis 
in they can be written as two coupled QLE, 



dt 



dt' [T„it ~ t')Q,{t') 



rr{t-t')Q2it')]^Bi{t) (1) 



and a similar equation where 1 and 2 are interchanged. 
Here, we introduced a distance d dependent damping ker- 
nel Tdit) = 70(e-"l*-''l -he-"l*+^l) , and bath operators 



where gk = {hgl/nikUJkm^y^^, and 6j.,6fe and a].,afc are 
creation and annihilation operators of a symmetric and 
antisymmetric bath mode k. Note that the operators 
Bi/2{t) evolve freely in time; the back-action of the two 
oscillators on the bath modes is solely contained in the 
memory terms in the QLE. The QLE also have a clear 
classical interpretation: the two oscillators are subjected 
to friction with a damping constant 7, they are coupled 
via a bath-mediated retarded interaction, and they are 
exposed to stochastic forces By^it). Without Vc the QLE 
also would exhibit terms proportional to Qi{t)Q2{t), cor- 
responding to an instantaneous, direct coupling of the 
two oscillators. In principle, the appearance of such term 
is possible because our model does not obey Lorentz 
invariance. Nevertheless, here we are interested in the 
bath-mediated coupling of the oscillators, and therefore 
eliminated the direct couplings by adding Vc to the sys- 
tem Hamiltonian. 

The formal solution of the QLE is simple, once they 
are written in the form 



Zy{t) + j^J^dt'C{t-t')y{t')^B{t) 



(2) 



where y ^ {QuQ2.Q1.Q2), B = (0,0,Bi,B2), and Z 
and C{t) are 4x4 matrices whose definitions become 
obvious by comparison of Eq. ([2|) with the original QLE. 
Then, the solution y(t) of Eq. ([2]) for initial y(0) and 
inhomogeneity B(i) is 



y{t) = g{t)y{0) 



dt' g{t - i')B(i') 



(3) 



where the Green's function Q{t) solves the homogeneous 
part of Eq. Its Laplace transform Q{s) = [s + Z + 
sC(s)]~^ can be calculated analytically. 

Correlations and entanglement in the two oscillator 
system can be studied on the basis of the oscillator's di- 
mensionless covariance matrix C, 

Cim = (y;ym + ymyi)p, = tr[(y;y,„ + y,„y;)Ps] , 

where ps is the joint state of the system oscillators. The 
vector y is obtained from y by multiplying the first and 
second entry with {mujQ/hy^^, and the third and forth 
entry with (m/fiwo)^^^^. Assuming that at time t = the 
total state factorizes in an initial oscillator state ps and 
a thermal state px of the bath, the temporal evolution of 
the covariance matrix follows with Eq. ([3]) to be 

c{t) = g{t)c{o)g{ty+ fdt' fdt"g{t-t')icit'-t")g{t-t"y . 

Jo Jo 

Here, C(0) is the covariance matrix of the initial oscilla- 
tor state Ps. The matrix JC{t) 2m{B{t)B{0y)p^/hwo 
contains the correlations of the bosonic fields B1/2. Its 
non- vanishing entries are IC^iit) = /C43(t), equal to 
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duj LU- 



■ coth 



2f 



■ cos ivt coswr , 



and two diagonal elements JC33{t) — JC^iit), which are 
given by the same expression, but with r = 0. 

We quantify the amount of entanglement of the two 
oscillators by the logarithmic negativity E. In case of 
a Gaussian state Ps of the oscillators, E can be conve- 
niently determined from the correlation matrix C as fol- 
lows: First, one applies a time-reversing operation [15[ on 
the second oscillator, according to which the covariance 
matrix transforms to Cim Cim — (— l)''"^*"■'C^m■ 
Then, the symplectic eigenvalues Ai, A2 of (7 yield the log- 
arithmic negativity as E = — X]j=i I0S2 -^j) UM- 
In this way the entanglement dynamics of the two os- 
cillators follows from the temporal evolution of the cor- 
relation matrix C, provided that the oscillator state is 
Gaussian for all times. This is the case when we restrict 
ourself to Gaussian initial states, since this property is 
conserved under the dynamics of the quadratic Hamilto- 
nian. Note that for Gaussian states a vanis hing logarith- 
mic negativity is equivalent to separability [l5|. 

Having outlined our model and the methods we have 
used, let us now present our results. Our main inter- 
est is the generation of entanglement from an initially 
separable oscillator state ps (0) via the coupling with the 
bosonic bath. Since there are no reasons for certain ini- 
tial separable states being preferred to other ones, here 
we present only results where initially the system oscilla- 
tors are in their ground state. The bath is assumed to be 
initially in a thermal state pr of temperature T. Thus, 
the total state is Gaussian and we can determine the log- 
arithmic negativity £' as a function of time as outlined 
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FIG. 1: Asymptotic entanglement of the system oscillators 
measured in logarithmic negativity i? as a function of distance 
r (in units of c/uq) for temperatures T = 0, 0.1, 0.2, and 
0.3 X hujo/kB (upper to lower curves), damping constant 7 = 
ujo, and cut-off frequency Q = lOuio- E drops to zero at a 
rather small critical distance do ^ c/fi, which for the above 
temperatures is proportional to (cf. inset). The fitted 
straight lines (dashed) have slopes a — 1.51, 1.18, 0.72, 0.25. 

above. In the following, we will mainly show numerical 
data demonstrating the characteristic dependence of the 
entanglement generation on distance r, cut-off frequency 
i7, damping constant 7, and temperature T. Generally, 
we measure distances in units of c/luq, frequencies in units 
of loq, and temperature in units of hioo/kB- 

First, we consider the entanglement of the two os- 
cillators at large times. For any finite distance r one 
finds Q{t) for i — » 00, meaning that the ini- 
tial oscillator state becomes irrelevant at large times. 
Hence, the asymptotic covariance matrix is Coo = 
J^dt'J^dt"g{t')IC{t' - t")g{t")'< . The time integrations 
together with the oscillating factors in IC{t' — t") rep- 
resent Laplace transformations, which eventually result 
in a single lo integral over terms containing the factor 
|0(ia;)p. The remaining integration over u can be eas- 
ily performed numerically. Fig. [T] shows the asymptotic 
logarithmic negativity i? as a function of the distance r 
between the oscillators. Clearly, the entanglement de- 
creases with distance and drops to zero at rather small 
critical distances do- The dependence of c?o on the in- 
verse cut-off frequency £7 for different temperatures can 
be seen in the inset of Fig. [TJ For Vl > ojq wc find the 
critical distance do to be inversely proportional to the 
cut-off frequency, dg ~ ac/fl, where a is a coefficient of 
order unity (at T ~ 0) that decreases with increasing 
temperature. The distance dg is rather insensitive to the 
actual value of the damping constant 7. For instance, 
the critical distance of do = 0.151 (in units of c/ujq) at 
7 = Wo and T = just changes to 0.12 or 0.17 when 
the damping is increased or lowered by a factor of 10, 
respectively. 

Now we consider how the logarithmic negativity de- 
velops in time. Determining the time dependent covari- 
ance matrix C{t) involves an inverse Laplace transfor- 
mation, which we performed numerically using Durbin's 
formula [l6| . Results for vanishing and three nonvanish- 
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FIG. 2: Logarithmic negativity E as function of time (in units 
of 1/ljo) for distances r = 0, 0.05, 0.1, and 0.15 x c/ujo (upper 
to lower curves) below the critical distance do, T — 0, Q — 
IOclio, and ^ = loq. Dashed lines represent asymptotic values. 
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FIG. 3: Short time behavior of E for r = 0.15, 0.2, and 
0.4 (in units of c/ujq, top to bottom). The initial peak of 
E{t) visible in Fig. [2] resolves in two peaks. The first peak 
is exponentially suppressed in 2rQ,/c. The second peak is 
delayed by approximately r/c. Its height decreases with r and 
vanishes for distances r > di. The inset shows di as function 
of the inverse cut-off frequency 1/fi (in units of I/ujq). 

ing distances r below the critical distance do are shown 
in Fig. [21 All curves show a characteristic peak at short 
times within which the logarithmic negativity reaches its 
maximum value Emax- After its decay the logarithmic 
negativity slowly recovers in an oscillatory manner to its 
asymptotic value, where the frequency of the oscillation is 
approximately wo/2. The oscillations decay rather slowly 
with time because the relative coordinate Qi — Q2 of the 
two oscillators is weakly damped for the small distances 
r under consideration. This behavior does not change 
much for distances slightly above dg. However, at larger 
distances r > O.ISc/cjq the logarithmic negativity does 
not recover at all but remains zero for all later times. 

Focussing on the short time behavior of E, the initial 
peak actually resolves in two peaks, as shown in Fig. [31 
The first peak appears immediately after switching on the 
interaction at times t less than r/c. Since bosons cannot 
have been exchanged between the two system oscillators 
within this time span, we attribute this peak to entan- 
glement that had been already present in the bath [l3]. 
Switching on the interaction might immediately transfer 
part of that entanglement to the oscillators. This behav- 
ior can be addressed by a short time expansion of C(t), 
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which, at zero temperature, eventually results in 

E{t) « 7^— (e-^ nt - ant [ntf + 0[mf] , (4) 

with an Vlt dependent a^t ~ 0.2937 — ^InVlt. ^From 
this expansion we find that for r > c/r2 width and height 
of the first peak are exponentially suppressed in the pa- 
rameter rf2/c by factors less than ~ exp(— rfi/c) and 
~ exp(— 2rf2/c), respectively. 

The second peak in the logarithmic negativity is de- 
layed by a little bit more than r/c, which suggests that 
it refers to entanglement due to exchange of bosons. Its 
height decreases monotonically with distance r and, in 
fact, reaches zero at a relatively small distance di that 
is constrained by the inverse cut-off frequency. At zero 
temperature and damping 7 = cjo numerical data show 
diin) < 6.0c/n (cf. inset of Fig.[S|. We expect that the 
actual value of the damping constant 7 has only minor in- 
fiuence on di (as like on the distance do), since numerical 
data as well as Eq. ^ show that in first approximation 
7 scales only the amplitude of E{t). 

We conclude that generally for distances r significantly 
larger than c/0 the logarithmic negativity E{t) reaches 
its total maximum Emax within an exponentially short 
time to < exp(— rri/c)/J7 and then vanishes for all times 
t > 2to. Moreover, at these distances the maximum value 
Emax is exponentially suppressed in 2rCl/c. 

To summarize, by analyzing the time-dependent loga- 
rithmic negativity of two oscillators coupled to a bosonic 
bath we found strong evidence that the entanglement 
mechanism under consideration is limited to rather small 
distances r of order of c/fi, i.e. to distances of order of the 
cut-off wavelength A. In practice, this length corresponds 
to the spatial extension of the microscopic systems to 
be entangled. At larger distances the maximum achiev- 
able logarithmic negativity is exponentially suppressed 
in, roughly, r/X. We believe that this behavior is charac- 
teristic for bath-mediated entanglement in general, since 
there seem to be no features of the investigated oscillator 
model which would it make special for entanglement. In 
fact, the general picture outlined here is fully supported 
by results that we obtained for an alternative two-spin- 
boson model [ist . Having said this, one may summarize 
our findings by stating that generally two objects can 
only be efficiently entangled via the interaction with a 
heat bath if they are in immediate vicinity of each other. 

It might appear puzzling that the environment quickly 
and strongly entangles with the two oscillators (which, 
after all, is the origin of the ubiquitous phenomenon of 
decoherence) , while the two oscillators for their own re- 
main essentially disentangled (if they are remote from 
each other). The reason behind this strongly asymmet- 
ric behavior is the large asymmetry in the (effective) 
Hilbert space dimensions of the participating systems: 



few oscillator states interact with a continuum of bath 
states. Assuming that the generic state of the joint sys- 
tem is well represented by a randomly chosen state of 
the joint system, it follows from [l^ that for dimensional 
reasons the bath is strongly entangled with each oscil- 
lator, while the system oscillators on their own remain 
separable. Thus, our analysis particularly demonstrates 
that under the actual dynamics - generated by a stan- 
dard bilinear system-bath interaction - a non-generic ini- 
tial state rapidly evolves to a generic one. Interestingly, 
the considered interaction fails to produce this effect if 
the distance r becomes less or of the order of the cut-off 
wavelength A, as evidenced in significant entanglement of 
the system oscillators in this case. 

We would like to thank C. Kiefer for valuable discus- 
sions. The work is supported by DFG grant TL2159 and 
by the Bonn-Cologne Graduate School. 
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